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Normal ground state of dense relativistic matter in a magnetic field 
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The properties of the ground state of relativistic matter in a magnetic field are examined within 
the framework of a Nambu-Jona-Lasinio model. The main emphasis of this study is the normal 
ground state, which is realized at sufficiently high temperatures and/or sufficiently large chemical 
potentials. In contrast to the vacuum state, which is characterized by the magnetic catalysis of chiral 
symmetry breaking, the normal state is accompanied by the dynamical generation of the chiral shift 
parameter A. In the chiral limit, the value of A determines a relative shift of the longitudinal 
momenta (along the direction of the magnetic field) in the dispersion relations of opposite chirality 
fermions. We argue that the chirality remains a good approximate quantum number even for massive 
fermions in the vicinity of the Fermi surface and, therefore, the chiral shift is expected to play an 
important role in many types of cold dense relativistic matter, relevant for applications in compact 
stars. The qualitative implications of the revealed structure of the normal ground state on the 
physics of protoneutron stars are discussed. A noticeable feature of the A parameter is that it is 
insensitive to temperature when T <C /io, where is the chemical potential, and increases with 
temperature for T > /io- The latter implies that the chiral shift parameter is also generated in the 
regime relevant for heavy ion collisions. 

PACS numbers: 12.39.Ki, 12.38.Mh, 21.65.Qr 



I. INTRODUCTION 

Dense relativistic matter in strong magnetic fields nat- 
urally exists in compact stars. For example, such type of 
matter is formed by the electron component of the nu- 
clear matter in the interior of neutron stars. At lower 
density, a relativistic electron plasma exists and plays an 
essential role in white dwarfs. In both cases, magnetic 
fields could be rather strong: they reach up to 10^ G in 
white dwarfs and up to 10^^ G in neutron stars [l|, Q. 
If quark stars exist in nature, the corresponding dense 
quark matter in the core will be a strongly coupled ver- 
sion of relativistic matter. Relativistic matter in a strong 
magnetic field is also created in heavy ion collisions (3| 
that can lead to the chiral magnetic effect d, • 

Many physical properties of the stellar matter under 
extreme conditions realized inside compact stars are un- 
derstood theoretically and could be tested to some extent 
through observational data. However, as was pointed out 
m Refs. the dense relativistic matter in a strong 

magnetic field may hold some new theoretical surprises. 
In particular, a topological contribution in the axial cur- 
rent at the lowest Landau level (LLL) was revealed in 
Ref. 0- More recently, it was shown in Ref. @ that the 
normal ground state of such matter is characterized by a 
chiral shift parameter A. The meaning of this parameter 
is clearest in the chiral limit: it determines a relative shift 
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of the longitudinal momenta in the dispersion relations 
of opposite chirality fermions, — > ± A, where the 
momentum k^ is directed along magnetic field. Taking 
into account that fermions in all Landau levels, includ- 
ing those around the Fermi surface, are affected by A, 
the corresponding matter may have unusual transport 
and/or emission properties. 

To further justify the motivation for this study, it is 
instructive to discuss the symmetry properties of the chi- 
ral shift parameter A. As we shall see below, it enters 
the effective Lagrangian density through the following 
quadratic term: A'07'^7^f/). Therefore, just like the ex- 
ternal magnetic field, the A term, being symmetric with 
respect to parity transformations V , breaks time rever- 
sal T and the rotational symmetry 50(3) down to 5*0(2) 
(i.e., the rotations about the axis set by the magnetic 
field). Also, since the A term is even under charge con- 
jugation C, it breaks CVT symmetry, which is also broken 
by the fermion density. We then conclude that the ab- 
sence of the chiral shift parameter is not protected by 
any symmetry, which, in turn, suggests that such a term 
should be dynamically generated even by perturbative 
dynamics. It is one of the purposes of this paper to shed 
light on this issue. 

The special role of the chiral shift parameter A will 
be discussed in detail below. Already here, however, we 
would like to point out that the quadratic part of the 
Lagrangian density A'07^7^?/; suggests a possible con- 
nection between the parameter A and the axial current 
along the direction of the magnetic field. Indeed, the 
parameter A enters the effective action as a Lagrange 
multiplier in front of the operator of the axial current 
j| = 'ip^'^^^ip. (This could be compared with the role 
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of the Dirac mass, which formally is the Lagrange mul- 
tiplier in front of the operator whose ground expectation 
value is the chiral condensate.) When the axial current 
is nonzero in the ground state, it should be generally ex- 
pected that A is also nonzero. Now, the axial current is 
known to be nonzero already in the system of noninter- 
acting fermions in an external magnetic field Q. Thus, 
we argued in Ref. Q that a nonzero A is an unavoidable 
consequence in interacting systems and, moreover, it is 
linear in the coupling constant to leading order [9|. In 
this paper we will confirm this suggestion. 

As was pointed out in Ref. 0, the structure of the 
topological axial current, induced only in the LLL, is in- 
timately connected with the axial anomaly [T^] ■ This fact 
is directly connected with the well known result that in 
a magnetic field the axial anomaly is also generated only 
in the LLL [l^. The important question is whether the 
form of the induced axial j| current coincides with the 
result in the theory of noninteracting fermions in a mag- 
netic field or whether it is affected by interactions (for 
related discussions, see Refs. 0, ajia, [H, [H, [ii). As 
has been recently shown in Ref. |2l|, while the dynam- 
ics responsible for the generation of the chiral shift A 
essentially modifies the form of this current, it does not 
affect the form of the axial anomaly. Moreover, while 
the topological contribution in the axial current is gen- 
erated in the infrared kinematic region (at the LLL), the 
contribution of A in this current is mostly generated in 
ultraviolet, which implies that higher Landau levels are 
important in that case. 

The main goal of this paper is to study in detail the 
dynamics responsible for the generation of the chiral shift 
parameter. This will be done at nonzero temperature 
and beyond the chiral limit in the Nambu-Jona-Lasinio 
(NJL) model. We also study some general and subtle 
features of the dynamics with the chiral shift parameter 
A. In particular it will be shown directly from the form 
of the gap equation in the NJL model that A necessarily 
exists in the normal phase in a magnetic field. Another 
property of A that is important for potential applications 
is that it is rather insensitive to temperature when T <C 
fiQ, where /io is the chemical potential, and increases with 
T when T > ^q. The first regime is appropriate for 
stellar matter, and the second one is realized in heavy 
ion collisions. 

Since the NJL model is nonrenormalizable, it is nec- 
essary to use a regularization with an ultraviolet cut- 
off. The important issue in such a model is how strongly 
the observables depend on the choice of a regularization 
scheme. In Ref. [9|, a gauge noninvariant regularization 
(with a cutoff in a sum over Landau levels) was used. In 
this paper, besides that regularization, we will also uti- 
lize the gauge invariant proper time regularization [2^ . 
It will be shown that the results in these two regulariza- 
tion schemes are qualitatively the same. 

The rest of this paper is organized as follows. In Sec. 
im the model of rclativistic matter in a magnetic field 
is introduced and the symmetry properties of various 



possible dynamical order parameters are overviewed. In 
Sec. Illli we derive the general form of the Schwinger- 
Dyson (gap) equation and discuss the approximations 
used in the analysis of the dynamics responsible for the 
chiral symmetry breaking and the generation of the chi- 
ral shift parameter. The analytical solutions to the gap 
equation at zero temperature are described in Sec. IIVI 
The numerical solutions of the gap equation are presented 
in Sec. [V] Both zero and nonzero temperature cases are 
analyzed in detail. The free energies of the corresponding 
solutions are calculated and the ground states for vari- 
ous sets of parameters are determined. In Sec. IVIl the 
induced axial current density in this model is calculated 
and analyzed. In Sec. I VIII we discuss our main results 
and their possible applications to the physics of compact 
stars and heavy ion collisions. Several Appendices at the 
end of the paper give many technical details and deriva- 
tions used in the main text. 



II. MODEL 

In this paper, in order to reveal the key elements of 
the dynamics responsible for the generation of the chi- 
ral shift parameter in the clearest way, we use the sim- 
plest Nambu-Jona-Lasinio model with one fermion fla- 
vor. Despite the obvious limitations, such a model with 
a short-range interaction is expected to provide a reason- 
able framework for revealing the qualitative features of 
the much more complicated dynamics in dense QED or 
QCD plasmas, where the long-range interactions are only 
partially screened. 

The Lagrangian density of the model reads 



C = i/j {iD^ + /io(5°) Yi' - moi^4> 

, Gint 



(1) 



where mo is the bare fermion mass and fxo is the chem- 
ical potential. By definition, 7^ = ij'^'y^j'^j^. The co- 
variant derivative Di, ~ ~ ieA^, includes the external 
gauge field A^. In the presence of a constant magnetic 
field pointing in the z-direction, the (3 -I- l)-dimensional 
Lorcntz symmetry in the model is explicitly broken down 
to the SO (2) symmetry of rotations around the z-axis in 
the presence of this magnetic field. Also, except parity 
V, all the discrete symmetries C, T, CV, CT, PT, and 
CVT are broken. 

In the chiral limit, mo = 0, this model possesses the 
chiral U{1)l x J7(l)fl symmetry. In the vacuum state 
(/io = 0), however, this chiral symmetry is known to 
be spontaneously broken at any G^it > because of 
the magnetic catalysis phenomenon [23 . [241 . (For lattice 
studies of this phenomenon, see Ref. [2^.) In essence, 
such spontaneous breaking results from the enhanced 
pairing dynamics of fermions and antifermions in the in- 
frared. The enhancement results from the nonvanishing 
density of states in the LLL that is subject to an effective 
dimensional reduction D — ^ D — 2. (This is somewhat 
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reminiscent of the pairing dynamics at the Fermi sur- 
face of a degenerate electron gas in the Bardeen-Cooper- 
Schrieffer theory of superconductivity.) At a sufficiently 
large value of the chemical potential, the chiral symme- 
try is expected to be restored. As we shall see below, this 
is indeed the case, but the corresponding normal ground 
state is characterized by a nonzero chiral shift parameter 
A. 



III. GAP EQUATION AT WEAK COUPLING 
A. Structure of the gap equation 

As follows from the structure of the Lagrangian density 
in Eq. ([l} , the tree level fermion propagator in coordinate 
space is determined by 



-,^3^,3 



Too 



(2) 



where u — (t, r), while tt^ = id'' + eA'', with k = 1, 2 , 
and TT^ = id'^ = —idz are the canonical momenta j2a |. 
(Note that the components of the conventional gradient 
V arc given by dk ^ —9'^ and the components of the 
vector potential A are identified with the contravariant 
components of the vector potential A''.) In the rest of 
this paper, we use the vector potential in the Landau 
gauge, A = {0,xB,0), where B is the strength of the 
external magnetic field pointing in the z-dircction. 

As for the structure of the full fermion propagator, it 
is given by the following ansatz: 

iG-'^{u,u') ^ + /^)7° - (7r_L • 7) - TT V 

+ ifi-/^-/^ + A-/^-/^-m]6'^{u-u'). (3) 

This propagator contains two new types of dynamical 
parameters that arc absent at tree level in Eq. p. 
and A. From its Dirac structure, it should be clear that // 
plays the role of an anomalous magnetic moment. As for 
A, it is the chiral shift parameter already mentioned in 
the Introduction. Note that in 2 -I- 1 dimensions (without 
z coordinate), A7'^7^ would be a mass term that is odd 
under time reversal. This mass is responsible for inducing 
the Chern-Simons term in the effective action for gauge 
fields |27| , and it plays an important role in the quantum 
Hall effect in graphene [2^, [23| . 

It should be emphasized that the Dirac mass and the 
chemical potential terms in the full propagator are deter- 
mined by TO and fi that may differ from their tree level 
counterparts, toq and fiQ. While toq is the bare fermion 
mass, TO has the physical meaning of a dynamical mass 
that, in general, depends on the density and tempera- 
ture of the matter, as well as on the strength of interac- 
tion. Concerning the chemical potentials, it is /xq that is 
the chemical potential in the thermodynamic sense. The 



value of /X, on the other hand, is an "effective" chemical 
potential that determines the quasiparticle dispersion re- 
lations in interacting theory. 

In order to determine the values of the parameters 
TO, fi, A and fl in the model at hand, we will use 
the Schwinger-Dyson (gap) equation for the full fermion 
propagator. As described in Appendix [Cl utilizing the 
approach based on the effective action for composite op- 
erators [13, HH , one can show that in the mean-field ap- 
proximation it takes the following form: 

G-\u, u') = S^\u, u') - iGint {G{u, u) ~ 7''^G(u, u)-f^ 
- tr[G(it, u)] + 7^ tr[7^G(u, u)] } 6'^{ 

(4) 

The diagrammatic form of the gap equation is shown in 
Fig. [T] While the first two terms in the curly brackets 
describe the exchange (Fock) interaction, the last two 
terms describe the direct (Hartree) interaction. 

This matrix equation is derived in the mean-field ap- 
proximation, which is reliable in the weakly coupled 
regime when the dimensionless coupling constant 



47^2 



(5) 



is small, g <^ 1. Here A is an ultraviolet cutoff, and the 
coupling g is defined in such a way that g^r = 1 , where gcr 
is the critical value for generating a fermion dynamical 
mass in the NJL model without magnetic field. 

Of course, weak coupling is completely adequate for 
the analysis of the electron gas in the interior of neutron 
stars. As for the stellar quark matter, such an approxi- 
mation may at best provide only a qualitative description 
of the dynamics responsible for the chiral asymmetry in 
the ground state. Regarding the use of the mean-field ap- 
proximation, there is no reason to doubt that it should 
capture the main features of the dynamics, especially in 
the weakly coupled limit. 





FIG. 1: Diagrammatic form of the gap equation in the 
Hartree-Fock (mean-field) approximation. 



As one can sec, the right hand side of the gap equation 
depends only on the full fermion propagator G{u, u') 
at u' = u. This fact greatly simplifies the analysis. Of 
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course, it is related to the fact that we use the local four- 
fermion interaction. This feature will be lost in more 
realistic models with long-range interactions. 

The main disadvantage of the local four-fermion inter- 
action is a nonrenormalizability of the model. Therefore, 
the model in Eq. (jlj should be viewed only as a low- 
energy effective model reliable at the energy scales below 
a certain cutoff energy A. One may try to associate the 
value of the cutoff with a certain physical scale, e.g., the 
Debye screening mass in dense matter, or another charac- 
teristic scale provided by nonperturbative dynamics. To 
keep it general, we assume that A is a free parameter in 
the analysis below. 



B. Structure of solutions of the gap equation 

In this subsection we consider the general structure 
of the solutions of the gap equation. In particular, it 
will be shown directly from the form of the gap equation 
that in the normal phase in a magnetic field a) ^ and /io 
are different, and b) the shift parameter A is necessarily 
nonzero. 

As shown in Appendix [Cl in the mean-field approxi- 
mation utilized here, the Dirac structure of gap equation 
^ does not allow solutions with a nontrivial p.. While 



having [1 = simplifies the analysis, we should empha- 
size that fl may well be nonvanishing in more refined 
approximations and in models with other types of inter- 
actions [2^, [2^, HH . At the same time, as one learns from 
a similar analysis in graphene, a nonzero fi should not 
change the main qualitative features of the phase with 
an induced A [2i,|2i|. 

The explicit expression for G{u,u) is calculated in 
Eq. (jA26p in Appendix The result reads 



G(it, u) 



E 



27r?2 ^ J (27r)2 



dudk^ JC-V-^ + JC+V+e{n - 1) 



(6) 

where ? = 1/ ■\/|eB| is the magnetic length, 9{n — 1) = 1 
for n > I and 6{n — 1) = for n < 0. We also use the 
following spin projectors: 



Pi = i (1 ± zs^tS^) , 



(7) 



and the shorthand notation s± = sign(e_B). The func- 
tions [/„ and /C,f are defined in Eqs. (|A24j) and (|A28j) . re- 
spectively. For reader's convenience, here they are quoted 
only for the case of vanishing /i, which is of the main in- 
terest: 



A2 



3^2 



2n\eB\ =F 2.s_lA (m + k^-f') 7°] 



(w + [if ~ 2n\eB\ ~ (^s^A - ^m^ -f (k?')'- 



{uj + [if - 2n\eB\ - ( .s_lA + ^rn^ + {k^) 



(8) 
(9) 



By making use of these expression, the nth Landau level contribution to the fermion propagator can be cast in the 
following form: 



±-D± 



Ur,. 



UJ + fi — is±j 



(^s±A - i/m^+TW 



LU + fi - isj_j 



I72 (^s^A+yG 



{k^y 



(w + [if - (sj_A - + (fc3)2^ _ 2n\eB\ {w + [if - (s_lA + + {k^f'^ - 2n\eB\ 



where 



m + k^j^ r o 



H± = - 1 ± S_L 



(fc3) 



3^2 



(11) 



are the projectors on the quasiparticle states, whose 
energies are given in terms of either the sum or the 
difference of s±A and -^/m^ -|- {k^f. The projectors 
take a particularly simple form in the massless limit: 
T~Lm=o ~ \ [1 ± sj_ sign(fc^)7^] . In this case, H^^g al- 
most coincide (up to the sign of the longitudinal momen- 
tum) with the chirality projectors Vf = ^ (l ± sj_7^), 
used in Ref. For each choice of the signs of eB and 
k^ , the chirality of the states that correspond to projec- 
tors H^j^Q arc summarized in TableUl In fact, the precise 



relation between the two sets of projectors reads 



^m=0 



1 =F sign(fc'^). 



l±sign(fc3) + 



(12) 



By making use of this relation and taking into account 
that \/m? + (fc3)2 _^ |;j3| .^yi^en m — > 0, it is straight- 
forward to check that the propagator in Eq. (jH) in the 
massless limit takes exactly the same form as in Ref. [§] , 
i.e., 



G{u, u) = G^V- + {G-V- + GlV+) , (13) 

n=l 
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TABLE I: Chirality of the eigenstates that correspond to pro- 
jectors H^^Q for each sign of the longitudinal momentum . 



and (p))) . 







=0 




=0 






fc^ > 


fc^ < 


fc^* > 


sign(e_B) > 


L 
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sign(eB) < 
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L 


L 


R 



where 



i 1 6517° f dujdk 



27r 



(27r)2 

w + /.J ± [fc3 - A sign(eS)] 



[uj + - 2n|eS| - [A:3 _ A sign(eS)]2 ^ 
w + /i =F [fc^ + A sign(eS)] 



(w + /i)2 - 2n|eB| - [fc3 + A sign(eB)] 



2 ' 5 



(14) 



The opposite chirality fermions, described by such a 
propagator, are characterized by a relative shift of the 
longitudinal momenta, k'^ — > fc^ ± sj^A, in their disper- 
sion relations. 

Unlike the higher Landau level terms in the propaga- 
tor, the LLL contribution is rather simple. 



UJ + fi- s±A + yjm^ + (fc3)2 



+ /i - s_lA - ^m2 + (fc3)2 



In the LLL, sj^ A is a part of the effective chemical poten- 
tial /X — s_lA, and the two terms in Eq. can be associ- 
ated with the antiparticle (negative energy) and particle 
(positive energy) contributions, respectively. In order to 
avoid a potential confusion, let us also mention that, as 
seen from Eq. ([TU)) . the connection between "H* and the 
particlc/antiparticle states is not preserved in the higher 
Landau levels. 

As follows from Eq. ((T0| and Eq. ([T5|) . the poles of the 
full fermion propagator are at 



wo = + s_lA ± v/m2 + (fc3)2, (16) 
for the lowest Landau level, and at 



A± V7 



3^2 



-2n|eB|, (17) 



for higher Landau levels, n > 1. Note that all four com- 
binations of signs are possible for the latter. 

The general form of the gap equation at nonzero tem- 
perature is derived in Appendix [Cl see Eqs. (jC4p . (|C5| 



m 
A 



n^o - GintS, 



(18) 
(19) 
(20) 



The functions A., B and V on the right hand side of these 
equations are determined by the full fermion propagator 
as follows: 

A = -tr[70G(^.,u)] ^ if), (21) 
B = -tr[G(u,M)] = {i^i;), (22) 
V = -tr[7VG(^x,u)] ^(jf), (23) 

and have the meaning of the fermion charge density, the 
chiral condensate and the axial current density, respec- 
tively. The formal integral representations of these three 
functions are presented in Eqs. (jA37|) . (|A38|) . and (|A40|) . 
Two of them, B and 2?, contain ultraviolet divergences. 
In the next section, we study various solutions to the gap 
equations by using two regularization schemes: a gauge 
noninvariant one, with cutoffs in momentum integration 
and the sum over the Landau levels |9|, and the gauge 
invariant proper-time regularization [22|. 

Let us now consider the zero temperature normal phase 
in the chiral limit, when m ~ rriQ = and (ij'ip) = 0. It is 
realized when the chemical potential /ig > T'ldyn/V^ (see 
Ref. and Sec. IIVI below), where rridyn is a dynamical 
fermion mass in a magnetic field at zero chemical poten- 
tial and zero temperature. Let us analyze Eqs. (|18|) and 
(|20p in perturbation theory in the dimensionlcss coupling 
constant g defined in Eq. ([S]). In the zero order approxi- 
mation, we have a theory of free fermions in a magnetic 
field. To this order, ^ = /xo and A = 0. However, even in 
this case the fermion density (j^) and the axial current 
density are nonzero. The former can be presented 
as a sum over the Landau levels: 

fio\eB\ sign(/xo)|eB| 



if) 



27r2 7r2 
x0(|/io|-v/2n|eB|), 



and the latter comes entirely from the LLL Q: 

-eB 



{jDo 



2n' 



-Mo ■ 



(24) 



(25) 



(The overall minus sign is due to our convention for the 
electric charge of the electron [1^.) Then, to the next 
order in the coupling constant, one finds from Eq. ([20)) 
that A oc Gint{ji)o 7^ 0. Thus, in the normal phase 
of this theory, there necessarily exists a shift parameter 
A. In essence, the latter is one of the main results of 
Ref. 0. Let us also emphasize that A is generated by 
perturbative dynamics, which is directly connected with 
the fact that the vanishing A is not protected by any 
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symmetry (recall that C = +1, V = +1, and T = —1 for 
the axial current density j|, and beside parity V, all the 
discrete symmetries are broken in model ([T])). 

This result was obtained for the case of zero tempera- 
ture. As will be shown below in Sec. I VII the chiral shift 
parameter is rather insensitive to the value of the tem- 
perature in the regime of cold dense matter appropriate 
for potential applications in stars. In the case of heavy 
ion collisions, as we shall see, a temperature larger than 
the chemical potential may play an important role in en- 
hancing the chiral shift parameter. 

As one can see from Eqs. (HH), and n - 

fiQ cx G'int(j'')o 7^ 0, which implies that ^ and fig are 
different. The origin of this difference can be traced to 
the Hartrce terms in the gap equation [see the last two 
terms in Eq. 

This finding seems to be robust in the NJL model with 
a local four-fermion interaction and a chemical potential, 
associated with a global charge, such as a baryon (or lep- 
ton) charge for example. When the conserved charge is 
related to a gauge symmetry, as in the case of the elec- 
tric charge, the situation may be different. In that case, 
a neutrality condition imposed by the Gauss law takes 
place [13 . The latter is necessary for providing the ther- 
modynamic equilibrium in a system. This is likely to 
result in fj.^^') = ^q'^'' when fj.^^') is the chemical potential 
for electric charge. Note that usually there are chemi- 
cal potentials of both types in dense relativistic matter. 
While being of importance for potential applications in 
principle, we expect that this fact will not change our 
main conclusion regarding the chiral shift parameter. 

In conclusion, let us briefly discuss the following issue. 
One may think that when the fermion mass m is zero, the 
term with the chiral shift A is unphysical: in this case, 
it could formally be removed by the gauge transforma- 
tion i/) — )• e^^'^^^tp, ■tp — >• ■06*^'''^^. The point, however, is 
that this transformation is singular (anomalous). It fol- 
lows from the two facts: (i) as was already pointed out 
above, in the LLL, s±A is a part of the chemical poten- 
tial (see Eq. ((15])), and (ii) this hap pens because the LLL 
dynamics is 1 + 1-dimensional [2J| . It is well known that 
in 1 + 1 dimensions this transformation, which formally 
varies the value of the chemical potential, is anomalous 
(for a recent thorough discussion of this transformation, 
see Ref. fH). 



IV. ANALYTICAL SOLUTIONS AT T = 

To set up a benchmark for the numerical results, it is 
instructive to start the analytical analysis of gap equa- 
tions p8)) . p^ . and (|20|) at zero temperature. We will 
use two regularization schemes: (i) the gauge noninvari- 
ant one, with a sharp momentum cutoff, |/c'^| < A, in the 
integrals over (which are always performed first) and a 
smooth cutoff in the sums over the Landau levels (which 
are performed last), and (ii) the gauge invariant proper- 
time regularization. It will be shown that the results in 
these two regularizations are qualitatively the same. 



1. Analytical solutions in the momentum cutoff 
regularization 

Let us start from the first regularization. The smooth- 
ing function in the sums over the Landau levels is taken 
in the following form: 



K{n) 



sinh(A/M) 



cosh ( {A/SA)^/n/n^] + cosh(A/(5A) 



(26) 



where the cutoff value ricut is determined by the number 
of the Landau levels below the energy scale set by A, 
i.e., TT-cut = [A^/2|ei?|] with the square brackets denoting 
the integer part. The width of the energy window in 
which the cutoff is smoothed is determined by the ratio 
A/(5A, and when the value of the latter goes to infinity, 
the function K{n) approaches the step function. (In the 
numerical calculations below we use A/ 6 A = 20.) 

Let us now show that there are two qualitatively differ- 
ent solutions, which were previously reported in Ref. 

Solution of Type I. The first solution type corresponds 
to m 7^ and A = in accordance with the magnetic 
catalysis scenario in the vacuum [23| ■ By substituting 
m 7^ and A = into the general expressions (jA37p . 
(|A38p . and (|A40[) . we derive the following expressions 
for the functions appearing on the right hand sides of 
Eqs. (HU-dini): 



7 



B 



V 



(tt/) 



n=l 



2(7r02 



2A |/i| + ^^2„^2 



E 



111 



A + + m2 + 2n|e_B| 
+ 2n|eB| 



In , = 1^1 - y/m^ + 2n|eB| 

+ 2n\eB\ ^ 

sign(^) 



2(7r/)^ 



■vV^^(H-H). 



(27) 



(28) 
(29) 



Note that the expression for V is proportional to the 
LLL contribution to the fermion density and, as a result, 
vanishes when < \m\. In this solution with 

A = is consistent with gap equation (PU]). Then, the 
other two gap equations reduce down to fi = fj.o and 



respectively. To leading order in the coupling constant, 
the solutions for fj, and A are straightforward. 



mo 



2gm 

00 

+ 2 ^ K{n) In 



1 2A 
In — 

A+ 



A = 



l + g/{Mf 



(36) 
(37) 



i2 +2n|eB| 



+ 2n|eB| 



, (30) 



where we utilized the smooth cutoff function (|26p in the 
sum over the Landau levels. 

Solution of Type II. In the chiral limit, in addition to 
the solution with a nonzero Dirac mass m, the gap equa- 
tion also allows a solution with m = and a nonzero chi- 
ral shift parameter A. To see this, we derive the functions 
that appear on the right hand sides of the gap equations 
for this special case: 



which are derived under the assumption that the chemi- 
cal potential fiQ is not large enough fo r the first Landau 
level to start filling up, i.e., < y^2\eB\. When the 
chemical potential becomes larger, the result for /.j will 
get corrections, but the expression for A in terms of /i 
will keep the same form. Note that with the function 
K(n) given in Eq. ([26|) . one finds that 



J2 n{^/2n\eB\,A) = aAV|eS| 



(38) 



where a = 0(1) 



A 
B 
V 



/.i - s_lA sign(/i) 



Nb 



2{ttIY 



J2 Vt^' - 2n\eBl (31) 

2. Analytical solutions in ttie proper-time regulanzation 



n=l 



= 0, 



(32) 



2{TTiy 



-A~2AY,^{n)], (33) 



n=l 



where Nb is the integer part of /{2\eB\). (Here it 
might be appropriate to note that the above result for T> 
remains unchanged also at nonzero temperatures!) The 
fact that now ;B = is in agreement with Eq. (|19p and the 
assumption m = mo = 0. The remaining two equations, 
(dH]) and (1201), reduce down to 



Solution of Type I. In the regime of magnetic catalysis, 
we have shown above that the dynamical mass parameter 
satisfies Eq. (|30p in the momentum cutoff rcgularization 
scheme. Now, let us show that this is consistent with the 
result obtained in the proper-time rcgularization. The 
expression for the vacuum part of function B in this reg- 
ularization is given in Eq. (|B9P . To the leading logarithm 
order, then, we derive the following gap equation for the 
mass parameter: 



and 
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{Aiy 



A - M - 2 sign(/i) - 2n\eB\^ 



mo 



mg 



In 



(AZ)2 7r(m0 



2 " 



(39) 



A 



{Aiy 



s_i_fi - A - 2A ^ K{n) 



(34) 



(35) 



As is easy to check, this gap equation is equivalent to 
Eq. (|30p to leading order. The corresponding gap equa- 
tion is also in agreement with the proper-time result in 
Ref. (23 | , which was given in the form 



m = mo + 



gm 



— e-^™'coth(|eB|s), (40) 

1/A2 S 



8 



where the proper-time cutoff sq = 1 / is convcntionahy 
given in terms of an uhraviolet energy scale A. By noting 
that coth(|eB|s) = 1 + 2 XI^^i e-^^l'^'^l and using the 
table integral 



ds 



1/A2 



r(o,^) 



In. 



+ 0(1), (41) 



we also sec that, up to higher order corrections in powers 
of 1/A^, the representation in Eq. (^0]) is equivalent to 
Eq. (1201) • 

As follows from Eq. (p9| . the value of the dynamical 
Dirac mass mdyn in the chiral limit, mo = 0, reads 



2J, 



dyn 



■ exp 



(42) 



Formally, this solution exists for |^ol < '^^dyn- As we 
will discuss below, however, it corresponds to the ground 
state only in a part of this range, |/j,o| ^ TOdyn/v^. 

Solution of Type II. Now let us consider the chiral limit 
and search for a solution with m = and a nonzero chiral 
shift parameter A using the proper-time representation. 
In this case, B = and the expressions for A and V 
are derived in Appendix [B) Function A is finite and, 
therefore, is given by the same expression as in Eq. ([5T|) . 
Function V contains ultraviolet divergences. The cor- 
responding regularized vacuum contribution is given in 
Eq. (jBlSp . By adding also the (finite) matter contribu- 
tion, derived in Eq. (|B20|) . we obtain 



V = 



-(A/A)=e,fi 



coth 



eB 
A2 



(43) 

where erfi(a;) = —ieri^ix) is the imaginary error function. 
By expanding the expression for 2? in inverse powers of 
A, we arrive at the following approximate result: 



2{t:IY 



(AO^ 



A 



(44) 



which is in agreement with the result in Eq. (j33p after 
making the identification ^(AZ)^ = 1 + 2^^]^K(n) ~ 
2a(AZ)^, where the parameter a is defined in Eq. (p8| . 
As we see, a = 1/4 in the proper-time regularization. 

The gap equation for /.t is insensitive to the ultraviolet 
dynamics and coincided with Eq. p4p . By making use of 
the approximation in Eq. (j44p , we arrive at the following 
equation for A: 



A 



(AO^ 



(AO' 



(45) 



which is equivalent to Eq. psp after the same identifi- 
cation of the regularization schemes is made. Also, the 
proper-time solution, 

Mo 



/i 



1 



{MY + h{My 



(46) 
(47) 



is equivalent to the solution in Eqs. ([55)1 and (P7)) . 



3. Free energy 

As should be clear from the above discussion, in the re- 
gion ImoI < "idyn, the two inequivalent solutions coexist. 
In order to decide which of them describes the ground 
state, one has to compare the corresponding free ener- 
gies. The general expression for the free energy density 
is derived in Appendix |DJ For the two cases of interest 
here, the corresponding results are given in Eqs. (|D21[) 
and ((D25|) . 



Myn 



2(27r/)^ 



(l + (mdy„0'ln|AZ|) (48) 



and 



i^A ^ ~ ,^ I 1 - 5- 



{2t:IY 



A2 



(49) 



respectively. In deriving the last expression, we used the 
approximate relations fjL c^i and A ~ gfiQeB/h?. By 
comparing the free energies in Eqs. (|48| and (|49|) . we 
see that the ground state with a nonzero A becomes fa- 
vorable when /.to ^ mdyn/'v/Z- This is analogous to the 
Clogston relation in superconductivity [ssj . 



V. NUMERICAL SOLUTIONS TO THE GAP 
EQUATION 

In order to solve numerically the set of gap equations 
(HH), (|T9| . and (|20|) . we have to regulate the divergences 
that appear in the integrals over the longitudinal mo- 
mentum and the sums over the Landau levels in the 
expressions for the chiral condensate B and the axial cur- 
rent density T>. In Sec. IIVI we used two regularizations: 

1) with a sharp momentum cutoff, < A, in the in- 
tegrals over fc"^ (which are always performed first) and 
a smooth cutoff in the sums over the Landau levels, and 

2) the proper-time regularizations. Because it was shown 
that at zero temperature the results in these two regular- 
izations arc qualitatively similar, we perform a detailed 
numerical analysis of the gap equations at arbitrary tem- 
perature by using the first regularization only, which is 
technically much simpler to implement. 

The form of the smoothing function K{n) in this reg- 
ularization is given in Eq. p6| . The width of the en- 
ergy window in which the cutoff is smoothed is deter- 
mined by the ratio A/(5A. When the value of this ra- 
tio goes to infinity, K{n) approaches a step function, 
^('T^cut ~ n)^ corresponding to the case of a sharp cut- 
off at ricut = [A^/2|ei3|]. We note, however, that taking 
a very sharp cutoff in the sums over the Landau lev- 
els may result in some unphysical discontinuities in the 
physical properties of the model as a function of the mag- 
netic field. This is because of the discontinuities in the 
dependence of the function ncut(|e^|), which defines the 
number of the dynamically accessible Landau levels. In 
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our numerical calculations, we choose a reasonably large 
value A/(5A = 20. 

In order to keep our model study as general as possible, 
we specify all energy/mass parameters in units of the 
cutoff parameter A. In the numerical calculations below, 
we use the following values of the coupling constant and 
the magnetic field: 



GintA^ 



= 0.25, 



\eB\ 



0.125A^ 



(50) 
(51) 



The coupling is rather weak to justify the approximations 
used in the analysis. In real dense or hot quark matter, 
the actual value of dimensionless coupling may be even 
stronger. In the degenerate electron gas in the interior of 
compact stars, on the other hand, it is still much weaker. 
Our purpose here, however, is to perform a qualitative 
analysis of the model and reveal the general features of 
the dynamics relevant for the generation of the chiral 
shift parameter. Therefore, our "optimal" choice of g is 
sufficiently weak to make the analysis reliable, while not 
too weak to avoid a very large hierarchy of the energy 
scales which would make the numerical analysis too diffi- 
cult. Similar reasoning applies to the choice of the mag- 
netic field in Eq. ([51]) . This is a sufficiently strong field 
that makes it easier to explore and understand the quali- 
tative features of the dynamics behind both the magnetic 
catalysis and the generation of the chiral shift parame- 
ter. In applications related to compact stars, the actual 
fields might be considerably weaker. However, this value 
may in fact be reasonable for applications in heavy ion 
collisions 0, |3l- 



A. Numerical solutions at = 

The gap equation is solved by multiple iterations of the 
gap equations. The convergence is checked by measuring 
the following error function: 



3 

E 



.-if 



\ ^-^ maxfx 

> z— 1 ^ 



2 2 
z,n ' i.n- 



(52) 



where Xi = /i, A,m for i = 1,2,3, respectively. (In the 
case when both Xi^n and Xi^n-i vanish, the correspond- 
ing ith contribution to e„ is left out.) When the value 
of e„ becomes less than 10~^ (at T 7^ 0) or 10^^ (at 
T = 0), the current set of A„ and m„ is accepted 
as an approximate solution to the gap equation. Usually, 
the convergence is achieved after several dozens of itera- 
tions. In some cases, even as few as five iterations suffices 
to reach the solution with the needed accuracy. This is 
often the case when we automatically sweep over a range 
of values of some parameter (e.g., the temperature or the 
chemical potential) and use the solution obtained at the 
previous value of the parameter as the starting guess to 
solve the equation for a new nearby value of the same 



parameter. However, even in this approach, the required 
number of iterations may sometimes be in the range of 
hundreds. This is usually the case when the dynami- 
cally generated A and m have a steep dependence on the 
model parameters, which is common, e.g., in the vicinity 
of phase transitions. 

In order to set up the reference point for the nonzero 
chemical potential calculations, let us start by presenting 
the results for the constituent fermion mass as a function 
of the bare mass mo at /^g = 0. The corresponding zero 
temperature dependence is shown by the black line in 
Fig. [21 As expected, the mass approaches the value of 
TOdyn in the chiral limit (mo — >■ 0). For the model pa- 
rameters used in this paper it reads: 



TOdyn 



7.1 X IQ-^A. 



(53) 



In the same figure, we also plotted the results for 
several nonzero values of temperature. These results 
show that the value of the dynamical mass in the chiral 
limit gradually vanishes with increasing the temperature. 
Within our numerical accuracy, the corresponding value 
of the critical temperature is consistent with the Bardeen- 
Cooper-Schrieffer theory relation, Tc ~ 0.57mdyn- We 
also note that the results for fi and A are trivial at all 
temperatures when /ig = 0. 




0.4 0.6 

ffJo/'Wdyn 



1.0 



FIG. 2: (Color online) The dependence of the constituent 
mass m on the bare mass mo at /xo = for several fixed values 
of temperature: T = (black), T = 3.75 x 10"* A = 0.53mdyn 
(blue), T = 5 X IQ-^A = 0.70mdyn (red), T = 10~^A = 
1.4l7Tidyn (green), r = 2 X lO^^A = 2.81mdyn (dark brown), 
and T = 4 X lO^^A = 5.62mdyn (light brown). The insert 
shows the details in the rectangular area around the origin. 



Numerical solutions at T = 



The solutions of A vs ^0 and m vs ^0, obtained by the 
iterating the set of gap equations, are shown in upper and 
lower left panels of Fig. [31 respectively. There the solu- 
tions for several different values of bare masses mo are 
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shown. It might be appropriate to note here that, in the 
vicinity of the phase transitions, we had obtained a pair 
of solutions for each fixed value of niQ . The two different 
solutions are obtained by sweeping over the same range 
of the chemical potentials fiQ in two different directions: 
(i) from left to right and (ii) from right to left. When 
such a pair of solution, forming a small hysteresis loop, is 
observed, a first order phase transition is expected some- 
where within the loop. To determine the location of such 
a phase transition, the comparison of free energy densi- 
ties for the corresponding pairs of solutions is required. 
The general expression for the free energy is derived in 
Appendix [D] By calculating the corresponding expres- 
sion for each of the solutions around the hysteresis loop, 
we could point the position of the actual phase transition. 
An example of such a calculation is presented in the lower 
right panel of Fig. |3l There the free energies of the pair 
of solutions in the case mo — are shown. The blue solid 
line correspond to the solution with m = and A ^ 0. 
The free energy shown by the red dashed line represents 
the other solution, in which the Dirac mass m is nonzero 
and A is zero at small /io < Wdyn- The free energies of the 
two solutions become equal at /xo,cr ~ 0.73mdyn- This is 
where the first order phase transition occurs. Note that 
the numerical value of /io.cr is within several percent of 
the analytical estimate mdyn 

Concerning the solution for fi vs /xq, the results are 
always such that /i « to within a few percent. There- 
fore, the corresponding plot would give little informa- 
tion. In order to get a deeper insight into the deviation 
of /i from fiQ, we find it instructive to plot the result 
for the difference /.jq — /i instead. Note that as follows 
from Eqs. (fT8|l and ((2T|) . the latter is proportional to the 
fermion charge density. The result is presented in the 
upper right panel of Fig. [3l We see that fig — ^ is always 
positive, meaning that the value of ^ is slightly smaller 
than fj,o. 

By comparing that graph for — fi with the depen- 
dence of the chiral shift parameter A on in the up- 
per left panel of the same figure, we observe that they 
have the same qualitative behaviors. In particular, A is 
nonzero in the ground state only if the fermion charge 
density is also nonzero there. In other words, the shift 
parameter is a manifestation of dynamics in a system 
with matter. Note that when the numerical values of the 
model parameters are used, the results for /x and A be- 
come ^ ~ /zo — X^Q and A ~ F/xq for mo = 0, where 
X w 0.0295 and Y « 0.0252. 



C. Numerical solutions at T 7^ and /io 7^ 

Let us now proceed with the numerical solution of the 
gap equation at nonzero temperature. At vanishing value 
of /j,o, several results for the constituent mass have al- 
ready been presented in Fig. [5) The other two parame- 
ters, ^ and A, were identically zero in that special case. 
Here we extend the solutions to nonzero values of /iq- The 



numerical results for A vs /io and m vs /io are presented 
in Fig. m Note that the dependence /x — /io vs /xo (not 
shown in that figure) is similar to that of A vs /io at all 
temperatures. As should be expected, temperature sup- 
presses the dynamical fermion mass (see the right panel 
of this figure). However, the situation is quite different 
for the chiral shift parameter. As one can sec in the left 
panel of the figure, A is rather insensitive to temperature 
when T ^ fiQ, and increases with T when T > /io- This 
property refiects the fact that higher temperature leads 
to higher matter density, which is apparently a more fa- 
vorable environment for generating the chiral shift A. 
While the first regime with T ^ /io is appropriate for 
stellar matter, the second one with T > /xq (actually, 
T 3> /io) is realized in heavy ion collisions. As wc discuss 
in Sec. IVIl] below, the generation of A may have impor- 
tant implications for both stellar matter and heavy ion 
collisions. 



VI. AXIAL CURRENT DENSITY 

It is instructive to calculate the ground state expecta- 
tion value of the axial current density. As we see from 
Eq. ([23|) . it coincides with the function V, 

= ~ tr[^' j'G{u,u)]^V. (54) 

In the case of the vanishing Dirac mass, m = 0, an ex- 
plicit expression for T) within the momentum cutoff and 
the proper-time regularization schemes were presented in 
Eqs. pS]) and (|44|) . respectively. Both expressions can be 
written in the same form: 

(j|)c.-^[/x-2a.^A(A0^], (55) 

where a is a dimensionless constant of order 1, deter- 
mined by the specific regularization scheme. When the 
proper time is used, we find from Eq. (j44[) that a = 1/4. 
In the case of the cutoff regularization, it is defined by 
Eq. (|38|) . Note that qualitatively the same result is also 
obtained in the point-splitting regularization [2ll |. 

The first term in the parenthesis in Eq. ([55)) is the same 
topological term that was derived in the free theory in 
Ref. [A, while the second term is an outcome of interac- 
tions [9[. It is interesting to note that, by making use 
of the gap equation (PH)) for A, the result for the axial 
current can be also rewritten in an alternative form: 

, 2A A A2 

Gint 27r^ g 

While this may not be very convenient in the free theory, 
in which both the coupling constant g and the chiral shift 
A vanish, and the cutoff is formally infinite, it is helpful 
to get a deeper insight in interacting theory. 

Formally, the results for the the axial current either 
in Eq. (j55p or in Eq. (j56p appear to be quadratically 
divergent when A — >■ 00. It should be noticed, however. 
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FIG. 3: (Color online) The zero temperature results for the chiral shift parameter A (upper left panel), the mass m (lower left 
panel), and the difference fio — fJ, (upper right panel) on the chemical potential fio for several fixed values of the bare mass: 
mo = (black solid line), mo — 10~*A — 0.14mdyn (blue long-dashed line), mo = 3 x 10"^ A — 0.42mdyn (red short-dashed 
line), mo = 5 x lO^^A = 0.70mdyn (green dash-dotted line). The lower right panel shows the free energies for the solution with 
a nonzero dynamical mass (red dashed line) and the solution with a chiral shift (blue solid line) in the chiral limit, mo = 0, at 
zero temperature. 




FIG. 4: (Color online) The nonzero temperature results for the chiral shift parameter A (left panel) and the mass m (right 
panel) as function of the chemical potential for several fixed values of the temperature, T = (black), T = 1.41mdyn (blue 
lines), T = 2.81mdyn (red lines), and T = 5.62mdyn (green lines), and two values of the bare mass, mo = 0.14mdyn (solid lines) 
and mo — 0.70mdyn (short-dashed lines). 



that the solution to the gap equation, see Eq. ([57)) in the 
case of cutoff regularization and Eq. (H71) in the case of 
proper time regularization, is inversely proportional to 
A^, i.e., A gi^LcB/A?. Taking this into account, we 
see that the axial current density is actually finite in the 



continuum limit A — > oo, 

, ^, -eB A^ . -eB n 

Os -^M + a— A- / (57) 

Before concluding this section, let us also note the fol- 
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lowing expression for fcrmion number density: 

(/) = -tv[j"G{u,u)] =A. (58) 

The explicit form of the function A is derived in Ap- 
pendix El The corresponding result is complicated and 
adds no new information when the solution to the gap 
equation is available. Indeed, the fermion number den- 
sity can be conveniently rewritten in a simpler form by 
making use of the gap equation (|18|) for fi, 



, Mq - M _ Mo ^ M 

Gint 



27r2 



9 



(59) 



This shows that the result for this density is proportional 
to Mo ~ M presented earlier. 

In the case of a strong magnetic field, when the LLL 
approximation is appropriate, we find that the chiral shift 
parameter (and, thus, the axial current density) and the 
fermion number density are proportional to each other. 
This is apparent in Fig. [31 The underlying reason for this 
proportionality is the same (up to a sign) LLL contribu- 
tion to both functions A and V. Moreover, this property 
seems to be at least approximately valid in a general case. 
In turn, this suggests that a fermion number density and 
the chiral shift parameter are two closely connected char- 
acteristics of the normal phase of magnetized relativistic 
matter. 



VII. DISCUSSION AND SUMMARY 

A. Fermi surface 

The immediate implication of a nonzero chiral shift pa- 
rameter in dense magnetized matter is the modification 
of the quasiparticle dispersion relations, see Eqs. (jl6p and 
P?)) . These relations can be used to define the "Fermi 
surface" in the space of the longitudinal momentum 
and the Landau index n. Note that the quantity 2n\eB\ 
plays the role analogous to the square of the transverse 
momentum A:^ = (fc^)^ -I- (fc^)^ in the absence of the mag- 
netic field. Following the standard philosophy, we define 
the Fermi surface as the hypersurface in the space of 
quantum numbers n and k^, which correspond to quasi- 
particles with zero energy, i.e.. 



: fc3 = ±v/(m^"7[A)2~^, 



(60) 



> : fc^ = ± J ( - 2n\eB\ ± sj_A) - m2(61) 



(In the last equation, all four combinations of signs 
are possible.) In order to better understand the na- 
ture of quasiparticles at the Fermi surface, described by 
Eqs. (|60)) and (|6T|) . we recall that there are two types of 
quasiparticles. The Dirac structures of their wave func- 
tions are obtained by applying the projection operators 
in Eq. ([TT|) . In relativistic dense matter {fi ^ m), the 



corresponding states at the Fermi surface can be approx- 
imately characterized by their chiralities. This follows 
from the fact that \k^\ ^ m for a large fraction of the 
Fermi surface in Eq. (j6ip . except for the limiting values 
of n around nmax = [M^/(2|ei?|)] . At such large values 
of the relative momentum, the projection operators in 
Eq. ([TT]) are very closely related to the chiral projectors. 
Indeed, for |fc^| ^ m, the relation between the two sets of 
projectors is approximately the same as in the massless 
case in Eq. ([T2|) . Taking this into account, it is possible to 
define quasiparticles at the Fermi surface, which arc pre- 
dominantly left-handed or right-handed. Without loss of 
generality, let us assume that sign(eB) > 0. Then, the 
Fermi surface for the predominantly left-handed particles 
is given by 



(62) 



n^O: k^ = +v/(/i- s_lA)2 -m2, 
n > : P = +\/ (^^^ ~ 2n\eB\ ~ s_lA)^ - m2,(63) 
k^ = (Vm^ - MeB\ + sj_Ay - m2,(64) 

and the Fermi surface for the predominantly right-handed 
particles is 



: k^ = -V(/i-s_LA)2-m2, 



(65) 



>0: k^ = -\ (y^fi^ -2n\eB\~ sj_A) - to2,(66) 



k^ = +\ i - 2n\eB\ + s^A) -m^.{67) 



In the massless case, of course, this correspondence be- 
comes exact. Then, we find that the Fermi surface for 
fermions of a given chirality is asymmetric in the direc- 
tion of the magnetic field. In Fig. |5l we show a schematic 
distribution of negatively charged fermions and take into 
account that the parameter sj_A has the same sign as 
the chemical potential, see Eqs. p7p or (|47p . (A similar 
distribution is also valid for positively charged fermions, 
but the left-handed and right-handed fermions will inter- 
change their roles.) For the fermions of a given chirality, 
the LLL and the higher Landau levels give opposite con- 
tributions to the overall asymmetry of the Fermi surface. 
For example, the left-handed electrons in the LLL oc- 
cupy only the states with positive longitudinal momenta 
(pointing in the magnetic field direction). The spins of 
the corresponding LLL electrons point against the mag- 
netic field direction. In the higher Landau levels, while 
the left-handed electrons can have both positive and neg- 
ative longitudinal momenta (as well as both spin pro- 
jections), there are more states with negative momenta 
occupied, see Fig. |S1 If there are many Landau levels oc- 
cupied, which is the case when fj, ^ -^/leBj, the relative 
contribution of the LLL to the whole Fermi surface is 
small, and the overall asymmetry is dominated by higher 
Landau levels. In the opposite regime of superstrong 
magnetic field (if it can be realized in compact stars at 
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FIG. 5: (Color online) A schematic distribution of (negatively 
charged) particles in the ground state of dense relativistic 
matter in a magnetic field (pointing in the positive z direc- 
tion). The values of the quantum number n (Landau levels) 
are shown along the horizontal axis, while the longitudinal 
momenta are shown along the vertical axis. The colored bars 
indicate the filled states of given chirality. 



When the original trapped neutrinos gradually diffuse 
through the bulk of a protoneutron star, they can build 
up an asymmetric momentum distribution as a result of 
their multiple elastic scattering on the nonisotropic state 
of left-handed fermions (electrons or quarks). This is in 
contrast to the common dynamics of diffusion through 
an isotropic hot matter, which leads to a very efficient 
thermal isotropization and, therefore, a wash out of any 
original nonisotropic distribution of neutrinos [stI . |38| . 

It appears also very helpful for the new pulsar kick 
mechanism that the chiral shift parameter is not much af- 
fected even by moderately high temperatures, 10 MeV ^ 
T ^ 50 MeV, present during the earliest stages of pro- 
toneutron stars [3^. Indeed, as our findings show, the 
value of A is primarily determined by the chemical po- 
tential and has a weak/nonessential temperature depen- 
dence when fi ^ T. In the stellar context, this ensures 
the feasibility of the proposed mechanism even at the ear- 
liest stages of the protoneutron stars, when there is suf- 
ficient amount of thermal energy to power the strongest 
(with V ^ 1000 km/s) pulsar kicks observed [s^. Alter- 
natively, the constraints of the energy conservation would 
make it hard, if not impossible, to explain any sizable 
pulsar kicks if the interior matter is cold (T ^ 1 MeV) . 

Let us also mention that the robustness of the chiral 
shift in hot magnetized matter may be useful to provide 
an additional neutrino push to facilitate successful su- 
pernova explosions as suggested in Ref . [i^l ■ The specific 
details of such a scenario are yet to be worked out. 



C. Heavy ion physics 



all), only the LLL is occupied and, therefore, the overall 
asymmetry of the Fermi surface will be reversed. In the 
intermediate regime of a few Landau levels occupied, one 
should expect a crossover from one regime to the other, 
where the asymmetry goes through zero. 

B. Effects in compact stars 

The asymmetry with respect to longitudinal momen- 
tum of the opposite chirality fermions in the ground 
state of dense magnetized matter, discussed in the pre- 
ceding subsection, may have important physical conse- 
quences. For example, the fact that only the left-handed 
fermions participate in the weak interactions means that 
the neutrinos will scatter asymmetrically off the matter, 
in which the chiral shift parameter is nonvanishing. 

By making use of this observation, a qualitatively new 
mechanism for the pulsar kicks |36| was proposed in 
Ref. 0. It can be realized in almost any type of rela- 
tivistic matter inside a protoneutron star (e.g., the elec- 
tron plasma of the nuclear/hadronic matter, or the quark 
and electron plasmas in the deconfined quark matter), in 
which a nonzero chiral shift parameter A develops. 



It is natural to ask whether the chiral shift parameter 
can have any interesting implications in the regime of rel- 
ativistic heavy ion collisions. As was recently discussed 
in the literature, hot relativistic matter in a magnetic 
field may have interesting properties even in the absence 
of the chiral shift parameter. The examples of the re- 
cently suggested phenomena, that appear to be closely 
related to the generation of the chiral shift, are the chi- 
ral mag netic effect 0, ; the chiral magnetic spiral 
[isl . ll4E6j , and the chiral magnetic wave [3l ■ 

As we find in this study, at high temperatures, i.e., in 
the regime relevant for relativistic heavy ion collisions, 
the chiral shift parameter is generated for any nonzero 
chemical potential. This is seen from the results pre- 
sented in Fig. m However, its role is not as obvious as 
in the case of stellar matter. At high temperatures, the 
Fermi surface and the low-energy excitations in its vicin- 
ity arc not very useful concepts any more. Instead, it 
is the axial current itself that is of interest. The chiral 
shift should induce a correction to the topological axial 
current (^5]) . As seen from Eq. ([55]) . the corresponding 
correction in the NJL model studied here is proportional 
to the chiral shift parameter A, multiplied by a factor 
(AZ)^, where A is the ultraviolet cutoff. Formally, the 
product of A and (AZ)^ is finite and is proportional to 
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the chemical potential. However, unlike the topological 
term, which is also proportional to the chemical poten- 
tial, the dynamical one contains an extra factor of the 
coupling constant. Therefore, only at relatively strong 
coupling, which can be provided by QCD interactions, 
the effect of the chiral shift parameter on the axial cur- 
rent can be substantial. 

Following the ideas similar to those that were used in 
the chiral magnetic effect would like to sug- 

gest that the axial current by itself can play an important 
role in hot matter produced by heavy ion collisions. It 
can lead to a modified version of the chiral magnetic ef- 
fect, which does not rely on the initial topological charge 
fluctuations. This can presumably be realized as follows. 
An initial axial current generates an excess of opposite 
chiral charges around the polar regions of the fireball. 
Then, these chiral charges trigger two "usual" chiral mag- 
netic effects with opposite directions of the vector cur- 
rents at the opposite poles. The inward flows of these 
electric currents will diffuse inside the fireball, while the 
outward flows will lead to a distinct observational signal: 
an excess of same sign charges going back-to-back. 

Concerning the regime of hot relativistic matter, let 
us also mention that it will be of interest to extend our 
analysis of magnetized relativistic matter to address the 
properties of collective modes, similar to those presented 
in Ref. [l^, by studying various current-current correla- 
tors. 



D. Renormalizability vs nonrenormalizability 

The present analysis was realized in the NJL model. It 
would be important to extend it to renormalizable field 
theories, especially, QED and QCD. In connection with 
that, we would like to note the following. The expression 
for the chiral shift parameter, A ~ gfieB/A?, obtained 
in the NJL model implies that both fermion density and 
magnetic field are necessary for the generation of A. This 
feature should also be valid in renormalizable theories. 



As for the cutoff A. it enters the results only because of 
the nonrenormalizability of the NJL model. 

Similar studies of chiral symmetry breaking in the vac- 
uum (/iQ = 0) QED and QCD in a magnetic field show 



that the cutoff scale A is replaced by \/\eB\ there 41 



Therefore, one might expect that in QED and QCD with 
both ji and B being nonzero, A will be replaced by a 
physical parameter, such as •^/|ei?|. This in turn suggests 
that a constant chiral shift parameter A will become a 
running quantity that depends on the longitudinal mo- 
mentum and the Landau level index n. 

Another important feature that one could expect in 
QCD in a magnetic field is a topological contribution in 
the baryon charge Q connected with collective massless 
fermion excitations in the phase with spontaneous chi- 
ral symmetry breaking. This feature could dramatically 
change the properties of that phase [TTI |. 

It is clear that dynamics in dense relativistic matter in 
a magnetic field is rich and sophisticated. In particular, 
one could expect surprises in studies of the phase diagram 
of QCD in a magnetic field 
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Appendix A: The full fermion propagator 
1. General result 

In this appendix we calculate the explicit form of the full fermion propagator. From the definition of the inverse 
propagator in Eq. it follows that [2Q| 



TO 



G{u, u') = i{u\ [{idt + m)7" - (^-L ■ 7) - t^^i'^ + iji-f'^^'^ + il^'^^^'^-f'^ 

= i{u\ [{idt + m)7*' - (7r± • 7) - 7r^7^ + ijl-y^^"^ - iA-f^j^-f'^ + ■ 
X I [{idt + m)7° - (tT-L ■ 7) - 7r^7^ + ifj."f^j^ + iA'y'^^^j^ - to] 

X [{idt + m)7" - (7r± • 7) - n^j^ + ifi^y^j'^ - iAj°-f^j'^ + to] ^ \u') 
= i{u\ [{idt + m)7° - i'^± ■ 7) - i'^7^ + nn^l^ - iA'j°"f^"f'^ + m] 

X [(i^t + /i + 1/2707172)2 - Tri - (7r^)2 - {m - tA-f°j^-/Y ' ieBj^j'^ - 2t-/'^-/'^{fl + Aj°)ttS^] \u'), 

(Al) 

where u = {t, r) and r = {x, y, z). Note that the canonical momenta are tt^ = id^ + eA^ (with fc = 1, 2) and tt'^ = — ic^a, 
and the Dirac structure 7'^7^ is rewritten in an equivalent form, 7'^7^ = z70'yi^^. 

In order to derive an expression for the propagator in a form of an expansion over the Landau levels, we follow the 
same approach as in Ref. [29[. Wc start by switching to the Fourier transform in time t and in z-coordinatc (i.e., the 
coordinate along the magnetic field), 

G(a;,fc3;r,r') = j dtdze''^^*-^''^-''''^'-'''^G{u,u') = i[W ~ {tt^^^ ■'y)]{r\{M - ir'i- ieB-i^-l'^Y^ \t'), (A2) 

where TTj^.r is the differential operator of the canonical momentum in the coordinate space spanned by vector r. The 
explicit structure of W and A^, which are matrices in Dirac space, follows directly from Eq. (|Aip . i.e.. 



W = {uj + ^)'^" + ifl-f^j^ - iAj°j^-f^ + m- k^-f^, (A3) 
M = {uj + ^i + iil-i°-i^-i^f -{m-iAf-i^-i^f -{k^f ~2i-i^-l'^{il + Af)k^-1^. (A4) 

By noting that all three operators A^, tt^ and ieB'y^'y^ inside the matrix element on the right hand side of Eq. (|A2[) 
commute, we proceed to build their common basis of eigenfunctions. We start from the operator tt^. As is well 
known, it has the eigenvalues (2n + l)|ei3| with n = 0, 1, 2, . . . . The corresponding normalized wave functions in the 
Landau gauge, A = (0, Bx), read 

V^„p(r) EE {r\np) = i_ J-^ ^^^ il + P^) e-i^^^+^'^^'e'^-™, (A5) 

where Hn{x) are the Hermitc polynomials [i^ and / = 1/ ^J\eB\ is the magnetic length. These wave functions satisfy 
the conditions of normalizability and completeness. 



Sv{np\v){r\n'p')= J d^rrnpirH'n'p'{r) = S,,,,,5{p - p'), (A6) 

00 

^ / dp{r\np){npW)=Y, / dp^np(r)Cp(r') = ^(r-r'), (A7) 

n n ^ 



n=0 



respectively. Then, by making use of the spectral expansion of the unit operator (jA7p . we can rewrite the matrix 
element on the right hand side of Eq. (|A2|) as follows: 



00 ^ 

{r\ {M - nl - ieB-f^-f'^y'^ \r') = ^ dp {r\np) [M - {2n + l)\eB\ - ieB-f^-f'^]~'^ {np\r') 

n 'J 

^i*(r,r') °° T (f\ 

.P-C/2V ^A8^ 

2ttP ^ M - (2n+ l)|eB| -16^7172' ^ 



Tl = 



16 



where ^ = (r — r')^/(2P) and $(r,r') is the Schwinger phase, whose exphcit form is given by 

r 

<I'(r,r') = -eJ dz^Mz) = -,s^i^±^2^^. (A9) 

r' 

This phase has a universal form for charged particles in a constant magnetic field. Its origin is related to the fact that 
the commutative group of translations is replaced by the noncommutative group of the so-called magnetic translations 
47| . In the derivation of Eq. (jASp , we calculated exactly the integral over the quantum number p by making use of 



formula 7.377 from Ref. l46j . 

oo 

e-^'i7,„(a; + y)ff„(a; + z)dx = 2"7ri/2m!z"-™L;;-"(-2y2), (AlO) 

which assumes m < n. By definition, are the generalized Laguerre polynomials, and L„ = L'^ [46j . 

By noticing that the matrix ieBj^j^ reduces down to its eigenvalues ±|ei3| in the subspaces defined by the spin 
projection operators, 

P± = i [l±Z7Vsign(eB)] , (All) 
the nth term in the sum in Eq. (jA8|) can be conveniently decomposed into the sum of two contributions, 



M - {2n + l)\eB\ - ieB-f^-f^ M - {2n + l)\eB\ + \eB\ M - {2n + l)\eB\ - \eB\ 



(A12) 



Here the ordering of the matrix factors 'P± and (M — 2n\eB\)^^ is of no importance because A4 commutes with both 
projectors. By substituting the last expression into Eq. (|A8[) and redefining the summation index n — > n — 1 in the 
second term, the result for the matrix element can be written in a compact form as 

n=0 ' ' 

where L_i = by definition. Finally, by noting that 

..-..(....^), 

the full propagator (jA2p takes the form of a product of the Schwinger phase factor and a translation invariant part, 
i.e., 

G(a;, r, r') = e'*(""'''')G(w, fc^; r - r'), (A16) 
where the translation invariant part of the propagator is 



G(w,fc^;r-r') 



M-2n\eB\ 

(A17) 

Note that the ordering of the matrix factors in this expression is very important because the expression in the square 
brackets does not commute with matrix A4. 

The Fourier transform of the translation invariant part of propagator (jA17[) . 

G(w,fc^k)= /d"re-'('^"")G(w,fc3;r), (A18) 
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can be evaluated by first performing the integration over the angle in the coordinate space. The integration is 
performed by making use of the following table integral: 



2ir 



j ^-^krcos{<P-<|,^)^^ ^ 27rJo(fcr), (A19) 


where Jo{x) is the Besscl function. Using also formula 7.421.1 from Ref. [iHl, one gets 
valid for i/ > and Re a > 0. The result is given by 

oo ^ 

G(a., k\ k) = ze-^-^'^ ^(-l)"i?„(a;, k\ k) ^ _ ^ (A21) 

where k = (fc^, k^) is the "transverse momentum", fc^ = |kp and the nth Landau level contribution is determined by 
Dr,{uj, k\ k) = 2W [P-L„ (2fc2p) - V+Ln-i (2fc2/2)] + 4(k . 7)L,\_i (2fc2/2) . (A22) 

It is assumed that j = 0. This could be compared with the standard Dirac propagator for massive fermions in a 
constant magnetic field in Ref. p^ . The last matrix factor in Eq. (jA2ip can be rewritten in a more convenient form 

as 

1 _ + - ijU7°7i72)2 _ (m + zA7°7i72)2 + 2i7i72(/i + A7°)fc373 - (fc3)2 _ 2n|eB| 

— , (A23) 



M - 2n\eB\ Un 
where 

Un = [(w + iif + A' - - A2 - (fc3)2 _ 2n|eS|]' - 4 {(Jl (w + + mAf + (k^f{A^ - . (A24) 

As follows from this representation, the fermion dispersion relations are determined by the zeros of J7„, which is a 
fourth order polynomial in in a general case. 

2. Fermion propagator in the coincidence limit, u' — > u 

The coordinate space representation of the propagator reads 

G(u, v!) = e'*'"-'--') / ^'■^^feY'^ e-'"(*"*')+''''<^--')+''--<'-^')G'(a;, k\ k). (A25) 
J (27r)4 

As seen from Eq. (U) in Sec. IIIIl only the full fermion propagator G'(m, u) in the coincidence limit, u' — >■ u, enters the 
mean-field gap equation. By making use of the results in the previous subsection, we obtain the following expression 
for the corresponding propagator: 

^, , /• dujdk'^d^k-. , 3 , , t ^ f dujdk^ JC-V- + JC+V+9{n - 1) , ^ ^. 

^^"''^) = j/ ^2^^("'^ '^)=2^gy -(2^^ it^ ^^''^ 

where 9{ti — 1) = 1 for n > 1 and 6{n — 1) = for n < 1. We took into account that $(r, r') = at r' = r, used the 
table integral 7.414 (8) in Ref. [i^, 

(valid when Re a > — 1 and Re s > 0) and introduced the following shorthand notation: 

IC^ = [{Lu + iJTs±A)"f° ±s^fl + m-k^-f^]^{uj + ^if +fl^ -A^ -{k^f -2n\eB_L\ 

T2s_l + + Am] 7° ± 2s±{il + Aj°)k^j^y (A28) 
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Note that the first factor in the square brackets originates from matrix W, see Eas. (jA3|) . (|A2ip and (jA22p . 

In the rest of this Appendix, we consider only the case with the vanishing anomalous magnetic moment (/i = 0), 
which is sufficient for the mean-field analysis in the NJL model. The expressions for /C„ and t/„ in this special case 
are quoted in Eqs. ([5]) and (O in the main text. Then, as follows from Eq. (|A26p . the propagator has the following 
general structure: 



Giu, u) - G^V- + J2 {Gn-P- + G+V+) 



(A29) 



By substituting /i = and replacing the integration over w by the Matsubara sum according to the prescription 



rn' — — oo 



with uj — > iujm' = iT^T{2m! + 1). we derive separate nth Landau level contributions, 

T 



G„ — — 



27r?2 
T 



-oo 27r Un 

^ m/ — — oo 



47r/2 



T 
4^ 



m — k y 



3„,3 \ lUJr, 



m' — — oo 



dk^ 



m — — oo 



where E^,^^ = 



sj^A ± -^Z m? + (fc3)2 j _^ 2n|ei3|. We use the following table sums: 

h 1 — np{b + ^) — np{b — ^) 



T E 

rn' — — ' 
oo 

^ E 



(a;,„' ~ ifi)^ + b'^ 2 

zw„i' + /.J np{b + ji) — np{b — ^) 



nv ——oo 



{lo^, - ifi)^ + 62 2 
where np^x) = l/(e^ + 1) is the Fcrmi-Dirac distribution function, and finally obtain 



1 



°° dk^ 



SttP ./_oo 27r 



1 



± 



87r/2 
1 

8^ 
1 

8^ 



27r 



7°± 



7°T 



— 70± 



27r 



°° ^1.3 



dy 

"27 



7°T 





(fc3)2 


m 




y^m? + 


(fc3)2 


m 




y^m'^ + 


(fc3)2 


m 




v/m2 + 


(fc3)2 



ym2 + (P)2 _ s_lA 



(A30) 



(A31) 



(A32) 
(A33) 



ym2 + (P)2 + s_lA 



l-np[ E^,^^ + ^L]-np[ E^, - fi 



l-np( , „ + - E.,^ - M 



(A34) 



The terms odd in k^ were dropped because they vanish after the integration over k^ is performed. 

In the case of the lowest Landau level (n = 0), the result can be rewritten in a more convenient form as 



G, 



1 



" 47r/2 
1 



dk^ 



+ 



2tt y/m^ + (fc3)2 
°° dP r 



l-np (^^/m'^ + (/e3)2 - sj^A + /i^ - n_F (^v^m2 + (fc3)2 + .s_lA - fij 



v/m2 + (fc3)2 - s±A + m) - np [^m^ + (k?Y + s_lA - /i 



(A35) 
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We notice that functions contain only two different Dirac structures: 7° and the unit matrix. Taking this into 
account, we conclude that the fermion propagator in the coincidence limit, v! — >■ u, is given in terms of just four 
independent Dirac structures: 

G(u, = YI'A + B + i-i-i^C + 7^7^!?] . (A36) 
The explicit expressions for functions A., B, C and V follow from the results in Eqs. (|A34p and (|A35|) . namely 



A 



{2t:IY 
1 



00 

00 



E 



B 



np 
dk^ 



(A37) 



(271^2 + (fc3)2 

- .00 ^^3 



E 



1 - np (Vm2 + (fc3)2 - s_L A + ^) - riF (Vm2 + (fc3)2 + a - 

1 - np (e-, +ti)-np - /i 



^m2 + (A;3)2 - s±A 



e: 



,n 

yjrn? + (fc3)2 + s_L A 

.n 



(A38) 



dk^ 



{2tiIY ^m2 + (fc3)2 



l-n_F (^v^77i2 + (A:3)2 - s_lA + ^ij - np (^^/mF+Jk^ + s_i_A - 



E 



dk^ 



i'2^^?^J-oo ^m^ + (fc3)2 



{2TTlf 



^ ^ n— 1 



E 



np {^-s/m? + (fc3)2 - A + A*) - '^f Vm2 + {k^Y + s_l A - ^ 

1 - np (e-s +h) - np - 



dk 



^Jm^ + (fc3)2 _ A 



^m2 + (fc3)2 + s_lA 



1 - ( Sj^3,„ + /i - ?iF I £^fc3.„ - M 



(A40) 



In order to clarify the physical meaning of these functions, it is useful to note their alternative definitions in terms of 
the following traces of the propagator: 



A = -tr[7°G(i.,zi)] ^(7^7V) = (/), 

B - -tT[G{u,u)] = {^y^), 

C ~ — tr [i7^72G(w, u)] = {ipi'f^'-f'^ip), 



(A41) 
(A42) 
(A43) 
(A44) 



Note that out of the four Dirac structures in the fermion propagator, only three of them A, B and V appear in the 
gap equation. Out of these latter, only two functions, i.e., B and V, contain ultraviolet divergences. In Appendix [B] 
below, we analyze functions A, B and V at zero temperature using the proper time regularization. 



Appendix B: Function A, B and D in proper time regularization 



1. Function ^ at T = 



As seen from Eq. (|A41[) . function A coincides with the baryon number density. As expected, the corresponding 
expression in Eq. (|A37[) has no divergences. Moreover, at T = 0, it can be given in a closed form in terms of elementary 
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functions, 



A\t=g 



sign(^ - s±A) 



sign(/i) 



Nb 



- s^AY - m^e - s^A| - |m|) 



2(7r02 



n=l 



^£y(V/^2_2n|ei?|-|A|) - sign (^A^^ - 2n|eS| - | A|) (^A^^ _ 2„|eS 



A| 



5] J VM2_2n|ei?| + |A 



V/i^ - 2n\eB\ + |A 



2(7r02 

where Nb is the integer part of ^'^ /{2\eB\). At m = 0, in particular, this becomes 



(Bl) 



A\T=0,m=0 



s±A sign(Aj) 



Nb 



(B2) 



2. Function S at T = 



As clear from Eq. (jA42p . function B coincides with the chiral condensate. We start by studying the divergent part 
of this function. As can be seen from Eq. (|A38p . all the divergences are independent not only of the temperature, but 
also of the chemical potential. Thus, at first we limit ourselves to the vacuum part (/i = 0) of the chiral condensate, 



dk 



3 



' V?7i2 + (fc3)2 - s_lA yjm? + (fc3)2 + s_lA' 



E 



n— 1 V '^k'^ ,n "^/c'^ ,n 



(B3) 



By reintroducing integration over frequency, this expression can be equivalently rewritten as follows: 



•'T=0 



/ ^m2 + (fc3)2 - .s_lA Jm? + {k^f + s^l\\ 

1^^ 



/"^ dk ^ ( Vm2 + k^ - s_lA ^ Vm2 + fc2 + s_lA 
2m fduj r „^ a;2+m2 + fc2- A2 + 27i|eS| 



/•°° dk ^ 



7 " J —Qo 



(27r02y TT [a;2 + (Vm2 + P _ g^i^y + 2n|eB|] p + (Vm2 + fc2 + s_lA)2 + 2n\eB\] 

(B4) 

where, by definition, a„ = 1 — ^iJj^. The divergent expression on the right hand side can be regularized by making 
use of the proper time method. In contrast to the cutoff rcgularization used in the main text, this regularization is 
explicitly gauge invariant. 

In Eq. (IB4[) . we introduce a proper time representation for each of the two factors in the denominator of the 
integrand and obtain the following result: 

oo 

2 iuJ^ + m2 + fc2 _ a2 + 2n|eB|)e-2"(^i+^^)l^^l 



cy ^ p J pOO pOO pOO 



Xg"(si + S2)(w +m" + fe" + A")g2siA(s2-,n)\/im2 + fe2 (^B5) 

Here it is convenient to introduce two new integration variables s and u instead of the proper times Si and S2, which 
are related as follows: 

si='^il-u), s2-^(l + u). (B6) 
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The integration over s runs from to oo, the integration over u runs from ~1 to 1, and the Jacobian associated with 
the change of variables is s/2. Then, summing over the Landau quantum number n and integrating over a; and u, we 
derive 



Rvac 2m 



OO /"OO 



dk 



ds sinh(2sAVm2 + P) 



— + + k'^ - coth(eBs) H ^ 

2s / sinh2(eBs) 



(B7) 

Note that in the last expression we introduced a finite proper-time cutoff at s = l/A" to regularize the ultraviolet 
divergences. It is easy to see that the divergent part at A — >■ oo is independent of A. Thus, after taking the limit 
A — > and integrating over k, we arrive at the following result: 



Kjvac 
"T=0 



ds 



- +m.^ ] coth(eB.s) H ^ 

s J sinh^(eBs) 



2(2^0' Ji/A 

- coth e-(™/'^)^ + 1^°^^^ ^ cotHeBs) 



(B8) 



where, after noting that ei?/ sinh^(e-Bs) = — ^ coth(ei?s), we were able to simplify the second term in the square 
brackets by integrating it by parts. Finally, we derive the result. 



Kjvac I 



m 



{2-kIY 



eB 



coth ( ^ ) e-™ + (AO 



2 -mVA^ _ 



(mO^Ei(mVA2) 



- Ill + 2 In r f 1 + ) + [mlf (l + In 2 - ln(,77i)') 

2A^TO m 



■In- 



1 



[ln(A0^ + hi2- 1] +0 



By making use of the following leading order approximation for the remaining proper time integration: 



(B9) 



dr 



cothr 



dr 



dr 



(r cothr - 1) + 0(e) 



a2Ei(ea2) + ln At + 21nr ( I + ^ ] + (l + In2-lna2) +0(e), (BIO) 
7ra^ \ 2 



where Ei(z) is the exponential integral function and we use the same identities for the generalized Riemann zeta 
fimctions as in Ref. [2J|. Note that Ei(z) ~ — Inz — 7^; + z + 0{z^) for z — > 0. 
The matter part of the function B at zero temperature reads: 



»2mattcr yy yyvac 



2{ttIY 



In 



\H - s_lA| + ^{^ - s_lA)2 - TV? 



^^(|/i-,s^A|-|m|)-ln 



|A| + VA^ 



-em-\m\) 



m 



E 



dk^ 



yjm? + (A;3)2 - s_lA 



E 



(Bll 



In the limit A — > 0, the result can be presented in an analytical form. 



Kpmattcr | 



m \fi\ + ^ - 



2{Triy 



/, /I I I i\ , \ri\ + J V? - m? — 2n\eB\ ^ ( . , , — ^ ^ . 

(ttO^ yjm^ + 2n|eS| V 

(B12) 
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3. Function T> at T = 

The vacuum part = 0) of the function D in Eq. (|A40p at zero temperature is given by 



7") vac 



{2TTiy 



v/m2 + (fc3)2 + s_lA -e[-^/m^^ + {k^Y - s^l^ 



E 



(27r02 



1 ^'m? + (A:3)2 - s_lA 1 + (fc3)2 _^ g^^^ 



E 



(27r02 y TT 



n=l 

oo oo 



V m2 + fc2 _ 5j_ A 



V m2 + A;2 + s±A 



u;2 + (Vm^ + F - s_L A)2 + 2neB uj^ + (Vm^ + P + s_l A)^ + 2nei3 



2A /"do; „ ^ 



+ fc2 _ ^2 _ ^2 _ 2neS 



(27rZ)2 7 TT 

Ya^ ^ dk dsi / ds2 (m2 + fc2 - _ a2 - 2ne5)e-2"('^i+^2)eB 



[^^ + (V"i2 + fc2 - s_lA)2 + 2neB] [uj^ + {sjm^ + P + s_lA)2 + 2neS] 

+ /"OO 



^g-(si+S2)CLj2+m2+fc2+A2)^2s^A(s2-si)Vm2+fc2 



(B13) 



We integrate over w, change the proper-time variables as in Eq. (jB6[) . and obtain 

/■+°° sinh(2sAVm2 + k'^)dk 



7^ vac 



1 

~2v^(27r02 

X e 



fc2 - i- - A2 ) eoth(eBs) - - 



2s 



eB 



sinh (eBs) 



(B14) 



In the hmit of to — > 0, in particular, this can be calculated analytically, 



T-)Vac I 

-'^T=oi™—>o 



2(27r 
A 



1 



coth(eBs) A^ + — coth(eBs) 



2.S 



ds 



2(27r/)2 Jo s 



°° ds 

— coth(ei3s) 



2(27r/)2 ^ 



■ erfi(ysA) coth(eBs) 



A 



e-(A/A)^ iierfi ^ coth ^ 



A 



eB\ AA 



A2 



1/A2 

A A^ 

(27r)2 3(27r)2 ' ^ V^^' ^2 



2A^ 



2(27r/)2-^ A 

where we first integrated over u and then integrated over k by using following table integrals: 

f + 00 



^^sinh(2afc)^_,fc2 



+ 00 



dk k sinh(2afc)e 



TT erfi 



„3/2 



'TT — crh 

ds 



where crfi(a) = — ierf(ia) is the imaginary error function, which has the following asymptotes: 

2a 



erfi(a) 



erfi(a) 



1 + — + — + — + (a^) 
3 10 42 ^ ' 



for a — > 0, 



1 3 15 

2^^ 4^^ 8^ 



oi-L 



, for a — > cx). 



(B15) 

(B16) 
(B17) 

(B18) 
(B19) 
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The matter part of the function I? at zero temperature reads: 



•T^mattcr 



'Dt=0 ~ 25^=0 



2s_L sign(^ - sj^A) r- --^ ^ 2sign(A) r— 



(27r02 

(2^ 



(27r02 



77126* (|A| - |m| 



^7772 + (A:3)2 _ s_lA 



ImI - 



fe3, 



V7772 + (F)2 +S±A 



E 



for 777 — >■ 0. 



(B20) 



Appendix C: Gap equation 

In order to derive the gap equation, it is convenient to use the formahsm of the effective action for composite 
operators [sH . In the mean- field approximation, the corresponding effective action T takes the following form: 



dt I d3r|(tr [G(7t, li)])^ - (tr [-i^'Giu^u)])' 



r(G) = -i Tr [LnG^^ + S-^G-l\ 

- tr[G(77,M)G(7i,7t)] +tr [7^0(1*, u)7^G(7t,u)]} , (CI) 

where u = it, r). The diagrammatic form of the equation is shown in Fig.[6l The trace, the logarithm, and the product 



o 







FIG. 6: Diagrammatic form of the expression for the effective action in the Hartree-Fock mean-field approximation. 

S~^G are taken in the functional sense. The gap equation is obtained by requiring that the full fermion propagator 
G corresponds to the variational extremum of the effective action, 5T /5G = 0, the explicit form of which reads 

G~^{u, u') = S-^{u, u') - iGint {G{u, u) - -f^G{u, 77)7^ - tr[G(7t, u)] + 7^ tr[75G(7i, u)]} 6*{u - u'). (C2) 

By making use of the ansatz ([3]) for the full fermion propagator, this gap equation takes the following form: 



{^I - Aio)7° + + iA^Vl^ - 777 + 7^0 = - 2Gint [7°^ + 7'7'2?] + GintS. (C3) 

Note that function C (anomalous magnetic moment) does not contribute to the right hand side of the equation. This 
is in accordance with the statement in the main text that the anomalous magnetic moment vanishes in the mean-field 
approximation used in our analysis of the NJL model. This matrix equation is equivalent to the following set of 
algebraic equations: 



1 „ , 



A 



I 



:Gint^5, 



777 = 7770 - GintS, 



(C4) 

(C5) 
(C6) 
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together with fl = 0. An alternative form of the same equations is given in terms of the baryon number density, the 
chiral condensate and the axial current density, 

M = /iO-^GintC?"), (C7) 

A = -^Gi„tOl>, (C8) 
m = mo - Gint (i/'^) • (C9) 



Appendix D: Free Energy Density 



In this Appendix, we derive the expression for the free energy density Q. It can be given in terms of the effective 
action F evaluated at the solution to the gap equation, n = —T/TV, where TV is a space-time volume. Taking into 
account the general form of the gap equation, the corresponding expression becomes 



1 



r = -i Tr 

By making use of the following Fourier transform of the Green's function G{u, u'): 



duj 



(Dl) 



(D2) 



we rewrite the effective action F as 

T = -iT 



duj ^ 
2^ 



-Tr 



lnG-^{uj) + - {S-\uj)G{uj) - 1) 



(D3) 



Then, by following the approach of Rcf. [29[ , we obtain the expression for the free energy density, 

duj f (Pk 



An J (27r)' 



rtr { [{lu - /io)7° + (k • 7) + ^'7'] G(w, k) + i} 



(D4) 



The propagator Gs(w, k^, k) is given in Eq. (|A2ip in Appendix 1X1 By making use of its explicit form and the table 
integral in Eq. (jA27p . we can calculate the following three integrals that contribute to the free energy density. 



dk^cfk 
(27r)3 

dk^d^]i 



(c^-/io)7°G,(c^,k) 
3 (k.7)G,(w,k) 



i2n) 
dk^d^k 



^ ' ri=0 



, (uj - Mo)7°W/ (an - ^7°) \V- + V+^n - 1)] 

, (a„ - 67°) 2n\eB\ [P_ + P+6>(n - 1)] 
C4 ' 



(D5) 
(D6) 



fcVGs(w,k) 



^ ' 11=0 



g f ^j^s [an fe7° - (c + d7°)fc^7^) [V- + -P+ejn - 1)] ^^^^ 



Un 



where a„, &, c and d are mutually commuting functions, defined by the following expressions: 

a„ = {uj + n)^ +fL^ -m^ - - {k^)^ -2n\eB\, 

c = — 2i7"'"7^/i, 
d = -2i"/^-/^A. 



(D8) 
(D9) 
(DIO) 
(Dll) 



[Note that the factor 6{n — 1) in Eq. (|D6|) is added for convenience; it is optional because the result is proportional 
to n and the n = term is vanishing anyway.] 
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By dropping an infinite divergent term which is independent of the physical parameters, from Eq. (jD4[) we derive 
the following expression for the free energy density: 



i f dujdk^ \ 



(fc3)2] [(,^-^o)(^ + M + s±A) + (fc3)2 



+2{uj - /io)(s_Lm - /i) + + Am] + 2s±fl{k^f{m - s±jl) + 2s±A{k^f{uj + fi + s±A) 
dujdk^ I 



2i -c-^ f diodk"^ r, 9 

-2{uj - tiQ)fl [fi{uj + ^i) + Am] + 2{k^f (A2 - /i^) | . 
Taking p, ~ 0, we have 

{uj - hq){uj + /X + S_lA) + (fc3)2 



9. 



2{2-KfP 



dujdk'^ 



2 A2 - _ 2n|eS|] [{lo - ^io)iLO + /i) + {k'^f + 2n\eB\\ 

(D12) 



(a; + - (s±A + y/m^Tlk^^ {cj + - (s±A - y/m?TJk^^ 



dujdk^ 



2(2^)3/2 _/ ^^2 + (P)2 



(g; - Mo)("^)^ + (fc^)^(^ + + s±A) _ (cj - f-io)im)^ + (fc3)2(a; + ^ + sj_A) 
[uj + ^if- (s±A + v/m2 + (P)2)' (w + /i)2- (s^A- Vm2 + (fc3)2)' 



(2 



(tj - ^o)(^^ + m) + (k^)^ + 2n|eB| (w - ^Jo)(u; + ^j) + (fc3)2 + 2n\eB\ 



E 



diodk^ 



(27r)3/2^^7 ^^2 + (fc3)2 



(fc3)2s_LA 



(fc3)2s_LA 



(D13) 



At T 7^ 0, the result for the free energy density reads 



T 



2(27r02 



m' — — OO 



+ iAto)(t^m' - - *sxA) - (fc^)^ +i/io)(w„' - - isxA) - (fc^)^ 



E 
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E/'^^^ E 
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(lj„j/ + ino){uim' - ifJ-) - {k^Y - 2n \eB\ + i/io)(Wm' - ifJ-) ~ (k^)^ - 2n\eB\ 



T 



E 



dk^ 



E 



m' — — OO 



(D14) 



where E^,^^ = ^m2 + (fc3)2± s_j_a and ^± ,„ 



2 + (fc3)2 ± sj^Aj + 2n\eB\ for n > 1. In order to calculate 

the Matsubara sums in this expression, we used the table sums in Eqs. (jA32p and (jA33p . as well as the following 
result: 



7n' — — oc 



(w,„' - i/i)2 + 52 



0-^ E 



m' ——00 



(m + ^J■o){^^nl' + m) - a2 - 52 

(Wm' - iAi)2 + 52 



^2 ^2 

[^^(6 + m) - n-pib - fi)] — — [1 - npib + - npib - ^J,)] . 



(D15) 
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Note that the vacuum subtraction in the first hne was performed. This is necessary in order to render the sum finite. 
In the limit T — 0, the above result reduces to 



Xt=o = (^(Ai + Mo)sign(/i)0(/i2-62) + 
After performing the Matsubara sums in Eq. (|D14[) . we obtain 



(Die) 



2(27r/)2 
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2(2^0' i?fe3,o 

/'CO 
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(/i + /io) sinh (/i/T) 



sinh 






cosh 




+ cosh (/i/r) 



cosh S+ „/T + cosh (/i/r) cosh £;-3.„/T + cosh (/i/T) 



E 



(fc^)^ + 27i|eS|+ i?+,J + 



1 g r°° dk^ 



{ky + 2n\eB\+[E-,\ - 



^ ' n=l^ fe'%n 

Note that we used the following identities: 

np{b — ii) — np{b + fi) 
1 — np{b — fi) — nir(5 + /i) 



(fc3)2g^A 
^777.2 + (fc3)2 

(fc3)2g^A 
^m2 + (fc3)2 



cosh ( S+ „/Tj + cosh (/i/T) 
sinh(i?^3.„/T) 



cosh LB-, „/T + cosh (/i/T) 



sinh(/i/T) 



cosh(Vr) +cosh(/i/T)' 

sinh(fe/T) 
cosh(6/T) + cosh(/i/T) ' 



(D17) 

(D18) 
(D19) 



1. Free energy density at m, fj, ^ and A = /i = 

In this case the free energy density (|D17p reduces to 



+ Mo 
i^^ir Jo 

1 



dk-^ 



sinh (/i/T) 



sinh {p-/T) 



(27r/)2 



cosh (Ea/T) + cosh (/i/T) cosh (En/T) + cosh (/i/T) 

1 [{k''y + Ei]smh{Eo/T) >22, 1 [(fc3)2 + 2n|eB| + ^2] sinh (£;„/T) 



Eo cosh (£;□ /T) + cosh (/i/T) 



cosh (£;„/T) + cosh (/i/T) 



where E^ = ^/77i2 + (A:3)2 _,_ 2n\eB\ for n > 0. 

In the limit of zero temperature, the above result for the final part of the free energy becomes 



1 f m 
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(D20) 



e{^? - - 27i|e-B|). 



(D21) 
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Note that, in the calculation, we subtracted an infinite constant term: 



2n\eB\ 



The sum over the Landau levels in the first line of Eq. ()D2ip still contains a logarithmic divergence, i.e. 



E- 



ln(AO^ 



In numerical calculations, we use of the smooth cutoff (|26p to regularize this expression. 



2. Free Energy Density at A,/i 7^ and m = /i = 

In this case the free energy density is 



+ Mo 
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(D24) 



where E^s ^ = -^Z (fc^ ± sj_A)2 + 2n\eB\. At T = 0, after doing the subtraction and the integration, we arrive at 



(27r02 



(2.0^ it 
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[1] P. M. Woods and C. Thompson, "Soft Gamma Repeaters 
and Anomalous X-ray Pulsars: Magnetar Candidates," 
in Compact Stellar X-ray Sources, edited by W. H. G. 
Lewin and M. van der Klis, (Cambridge University Press, 
Cambridge, 2006) pp. 547-586 astro-ph/0406133 . 

[2] S. Mereghetti, Astron. Astrophys. Rev. 15, 225 (2008). 

[3] V. Skokov, A. Illarionov and V. Toneev, Int. J. Mod. 
Phys. A 24, 5925 (2009). 

[4] D. Kharzeev, Phys. Lett. B 633, 260 (2006); D. Kharzeev 
and A. Zhitnitsky Nucl. Phys. A 797, 67 (2007); 
D. E. Kharzeev, L. D. McLerran, and H. J. War- 
ringa, Nucl. Phys. A 803, 227 (2008); K. Fukushima, 



D. E. Kharzeev, and H. J. Warringa, Phys. Rev. D 78, 
074033 (2008); Seung-il Nam, Phys. Rev. D 82, 045017 
(2010). 

[5] P. V. Buividovich, M. N. Chernodub, E. V. Luschevskaya 
and M. I. Polikarpov, Phys. Rev. D 80, 054503 (2009). 

[6] D. Ebert, K. G. Klimenko, M. A. Vdovichenko, and 
A. S. Vshivtsev, Phys. Rev. D 61, 025005 (1999); 

E. J. Ferrer, V. de la Incera, and C. Manuel, Phys. Rev. 
Lett. 95, 152002 (2005); E. J. Ferrer and V. de la In- 
cera, Phys. Rev. D 76, 045011 (2007); K. Fukushima and 
H. J. Warringa, Phys. Rev. Lett. 100, 032007 (2008); 
J. L. Noronha and I. A. Shovkovy, Phys. Rev. D 76, 



28 



105030 (2007). 

M. A. Metlitski and A. R. Zhitnitsky, Phys. Rev. D 72, 
045011 (2005). 

D. T. Son and M. A. Stephanov, Phys. Rev. D 77, 014021 
(2008). 

E. V. Gorbar, V. A. Miransky, and I. A. Shovkovy, Phys. 
Rev. C 80, 032801(R) (2009). 

A. Rebhan, A. Schmitt and S. A. Strieker, J. High Energy 
Phys. 01 (2010) 026. 

F. Preis, A. Rebhan and A. Schmitt, J. High Energy 
Phys. 03 (2011) 033. 

K. Fukushima and M. Ruggieri, Phys. Rev. D 82, 054001 
(2010). 

G. Basar, G. V. Dunne and D. E. Kharzeev, Phys. Rev. 
Lett. 104, 232301 (2010). 

K. Y. Kim, B. Sahoo and H. U. Yee, JHEP 10, 005 

(2010). 

D. E. Kharzeev and H. U. Yee. larXiv:10 12.6026 [hep-th] . 
I. E. Frolov, V. C. Zhukovsky and K. G. Klimenko, Phys. 
Rev. D 82, 076002 (2010). 

S. L. Adler, Phys. Rev. 177, 2426 (1969); J. S. Bell and 

R. Jackiw, Nuovo Cimento A 60, 47 (1969). 

J. Ambjorn, J. Greensite and C. Peterson, Nucl. Phys. 

B 221, 381 (1983); N. Sadooghi and A. J. Salim, Phys. 

Rev. D 74, 085032 (2006). 

V. A. Ruba kov, arX iv: 1005. 1888 [hep-ph]. 

D. K. Hong. Im^X^I OIO.3923 [hep-th]. 

E. V. Gorbar, V. A. Miransky, and I. A. Shovkovy, Phys. 
Lett. B 695, 354 (2011). 

J. S. Schwinger, Phys. Rev. 82, 664 (1951). 
V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Phys. 
Rev. Lett. 73, 3499 (1994); Phys. Rev. D 52, 4718 (1995). 
V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Phys. 
Lett. B 349, 477 (1995); Nucl. Phys. B 462, 249 (1996). 
P. V. Buividovich, M. N. Chernodub, E. V. Luschevskaya 
and M. I. Polikarpov, Phys. Lett. B 682, 484 (2010). 
Here our convention is such that the electric charge of 
the electron is — e where e > 0. 

A. J. Niemi and G. W. Semenoff, Phys. Rev. Lett. 51, 
2077 (1983); A. N. Redlich, Phys. Rev. Lett. 52, 18 
(1984); Phys. Rev. D 29, 2366 (1984). 
E. V. Gorbar, V.P. Gusynin, and V. A. Miransky, Low 
Temp. Phys. 34, 790 (2008). 

E. V. Gorbar, V. P. Gusynin, V. A. Miransky, and 

I. A. Shovkovy, Phys. Rev. B 78, 085437 (2008). 

J.M. Luttinger and J.D. Ward, Phys. Rev. 118, 1417 



(1960); G. Baym and L.P. Kadanoff, ibid 124, 287 (1961); 

G. Baym, ibid 127, 1391 (1962). 
[31] J.M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. 

D 10, 2428 (1974). 
[32] E. J. Ferrer and V. de la Incera, Phys. Rev. Lett. 102, 

050402 (2009). 

[33] J. L Kapusta, Phys. Rev. D 24, 426 (1981); 

V. P. Gusynin, V. A. Miransky and L A. Shovkovy, Phys. 

Lett. B 581, 82 (2004). 
[34] T. Kojo, Y. Hidaka, L. McLerran and R. D. Pisarski, 

Nucl. Phys. A 843, 37 (2010), and references therein. 
[35] A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962). 
[36] A. G. Lyne and D. R. Lorimer, Nature 369, 127 (1994); 

J. M. Cordes and D. F. Chernoff, Astrophys. J. 505, 315 

(1998); B. M. S. Hansen and E. S. Phinney, Mon. Not. 

Roy. Astron. Soc. 291, 569 (1997); C. Fryer, A. Burrows, 

and W. Benz, Astrophys. J. 496, 333 (1998); Z. Arzou- 

manian, D. F. Chernoffs, and J. M. Cordes, Astrophys. 

J. 568, 289 (2002); S. Chatterjee et al., Astrophys. J. 

630, L61 (2005). 
[37] A. Kusenko, G. Segre, and A. Vilenkin, Phys. Lett. B 

437, 359 (1998). 
[38] L Sagert and J. Schaffner-Bielich, J. Phys. G 35, 

014062 (2008); Astron. Astrophys. 489, 281 (2008); 

arXiv:astro-ph/0612776 
[39] D. Page and S. Reddy, Ann. Rev. Nucl. Part. Sci. 56, 

327 (2006). 

[40] C. L. Fryer and A. Kusenko, Astrophys. J. Suppl. 163, 
335 (2006). 

[41] V. P. Gusynin, V. A. Miransky and L A. Shovkovy, Phys. 
Rev. D 52, 4747 (1995); Nucl. Phys. B 563, 361 (1999); 
V. A. Miransky and I. A. Shovkovy, Phys. Rev. D 66, 
045006 (2002). 

[42] A. J. Mizher, M. N. Chernodub and E. S. Fraga, Phys. 

Rev. D 82, 105016 (2010). 
[43] M. D'Eha, S. Mukherjee and F. Sanfilippo, Phys. Rev. D 

82, 051501 (2010). 
[44] R. Gatto and M. Ruggieri, Phys. Rev. D 82, 054027 

(2010). 

[45] S. Fayazbakhsh and N. Sadooghi. [arXiv:1009.6125 ^ [hep- 

ph]. 

[46] L S. Gradshtein and L M. Ryzhik, Tables of Integrals, 

Series, and Products (Academic Press, Orlando, 1980). 
[47] J. Zak, Phys. Rev. 134, A1602 (1964). 



